546 


The Transmission of Electric Waves round the Earth. 

By G. N. Watson, So.D., I3.Sc, Professor of Mathematics in the 

University. of Birmingham. 

(Communicated by Prof. J. W. Nicholson, F.R.S. Received January 13, 1919.) 

On the hypothesis that the Earth consists of an imperfectly conducting 
sphere surrounded by infinite homogeneous dielectric, I have recently* 
obtained a complete solution (in a form adapted for numerical computation) 
of the problem of determining the effect at a distant point of the Earth's 
surface due to a Hertzian oscillator emitting waves of a definite frequency. 
Previous investigators had obtained approximations*)- (some of which were 
incorrect) to the dominant terms of the series which represents the effect 
due to the Earth, but the earlier approximations cease to be valid in the 
neighbourhood of the antipodes of the transmitter. On this hypothesis the 
absolute value of the Hertzian function (with the time-factor suppressed) is 
roughly proportional to (sin 6)~* exp ( — 23*94 V^#), where X is the wave¬ 
length measured in kilometres, and 6 is the angular distance from the 
transmitter. When 6 is nearly equal to tt, the factor (sin 0)~~% has to be 
suppressed. 

This formula does not agree with results obtained experimentally. 
The numerical factor 23*94 is much too large, so that, as 6 increases, the 
magnetic force decays much less rapidly than the theory indicates; and it 
has also been suggested on experimental grounds^ that the actual state of 
affairs is represented much more closely when the factor X~% is replaced by 
the factor 

To secure a considerable modification in the magnitude of the factor 23*94, 
a suggestion has been put forward by various physicists to the effect that 
the attenuated gases in the upper atmosphere have an appreciable con¬ 
ductivity, and so they tend to reflect electric waves. This hypothesis was 
first made in 1900 by Heaviside, and it has been developed more recently-by 
Eecles.§ 

* 4 Roy. Soc. Proc,’ A, vol. 95, pp. 83-99 (1918). 

t References to the more important work on the subject are to be found in my 
previous paper. A complete bibliography is given by Love, 4 Phil. Trans,’ A, vol. 215, 
p. 106 (1915). 

J See e.g.y Austin, ‘Washington Bureau of Standards Bulletin’ (October, 1911), who 
suggests 9*6 X~i in place of 23*94 A“i 

§ Eccles has worked out some of the mathematical consequences of the theory, 4 Roy- 
Soc. Proc,’ A, vol. 87, pp. 79-99 (1912). 
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It therefore seems desirable to investigate the consequences of the 
assumption that the Earth is surrounded by a concentric conducting layer 
at a considerable height. The height which has been suggested by Eccles 
is from 100 kilom. to 150 kilom. above the surface of the Earth. This com- 
paratively sharply defined region of the atmosphere which plays the pre¬ 
dominant part in signalling by night is called “ Heaviside's layer ” by Eccles. 
Signals during the day are affected by partial ionisation (due to sunlight) of 
lower regions of the atmosphere. 

A consequence of this theory is that the disturbance produced by 
the oscillator is confined to the neighbourhood of the Earth; this 
concentration obviously tends to increase the magnetic force at points 
on the Earth's surface. And it will, in fact, be shown that, if the 
Earth be regarded as a perfect conductor surrounded by a concentric 
perfect reflector, the expression for the Hertzian function consists of a finite 
number of oscillatory terms combined with an infinite series of negative 
exponentials. 

If the finiteness of the conductivity of the Earth and reflector is taken into 
account, the character of the finite series is modified; and it appears that 
the Hertzian function constructed on the reflector theory agrees quite well 
with the form suggested by experimenters. In particular, it seems impossible 
to account for the factor X~* on any theory of pure diffraction; and the 
presence of convergent waves is necessary to produce it. 

The results of this paper therefore tend to confirm the existence of the 
reflecting layer. A consequence of its presence is that it places grave 
obstacles in the way of communications with Mars or Venus, if the 
desirability of communicating with those planets should ever arise. 

2. In this paper the consequences of the ionisation theory are investigated 
by supposing that the effect produced by the regions of the atmosphere which 
tend to reflect the waves may be represented by taking the Earth to be 
surrounded by a concentric reflector of radius c ; it is at first supposed that 
the reflector is perfect, so that the tangential component of electric force is 
zero, and subsequently the theory is modified by supposing the reflector to 
have a finite conductivity; this procedure removes the grave analytical 
difficulties which would arise if the dielectric surrounding the Earth were 
supposed to have a conductivity which varied continuously. With the 
exception of the novelty involved in the introduction of c, the notation will be 
the same as in my previous paper. 

The radius of the Earth is denoted by a, the distance of the transmitter 
from the centre of the Earth is 5, and of course a^h<c. 

The transmitter, which is taken to be on the axis of harmonics, emits waves 
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of period 2 t t/co, and the components* of the electric force and of the magnetic 
force are (E r , E 0 , E^) and (H r , U d , H^); there are two constants! /3, y of the 
medium transmitting the waves such that the electric force and the magnetic 
force are determined in terms of the Hertzian function H by the equations 



H r = 0,. H 0 = O, 


1 3 2 (rXI) 

rb dr 80 9 

jr _ ft an 

a<f) “ b dd ’ 


e *=°> 


where H satisfies the wave equation 

(v 2 +/^)n = o, 

in which h 2 = — fiy. 

If the Hertzian function were unaffected by the presence of the Earth and 
of the reflector, its value would be Ii 0 = e~ tJcii fR, where E is the distance 
from the transmitter; and H 0 can be expanded in the forms 

n QO 

Hi = — 2 (2n +1 ) r n (kb)'f, l (kr) P„<», (r<5), 

/CTO n - o 

a 00 

H° - 2 (2rc + l)&(ir)fnWP.W, G>G, 

where and yfr n are expressible in terms of Bessel functions by the equations 

'xf/'n {x) — {^TTX)^ J n + | G), G) “ (| Hjj + | (®), 


the last function being the second of the two Hankel-Melsen functions, some¬ 
times called Bessel functions of the third kind. 

We now take into account the presence of the Earth and of the reflector: 
the values! of /3, 7 , h in the interior of the Earth are denoted by the symbols 
fii, 7 i, h) and the disturbance produced in the Hertzian function in the 
atmosphere is n^, while the disturbed function in the interior of the Earth is 
■Hi The appropriate series for H^ and H* are 


% 00 

rr; ^ (2%-f- 1) (J^t) -f- Cyfyn (G\)} 

rvTO n — o 

n CO 

n< =- y —7 2 (2 n +1) b n yfr n {hr) V n {fi\ 

ikro n _ o 

where the coefficients a n , b ny and c n are constants. 

* The time factor is supposed to be suppressed throughout. 

t These constants depend also on the frequency of the waves ; their precise values are 

Q cr + ie <o iixoa 

* = -e-’ 7 c’ 


where o- is the conductivity, e is the dielectric constant, fx is the permeability, and C is 
the velocity of light; the symbol /x is used elsewhere to denote cos d, where 6 is the 
angular distance from the transmitter, but this will cause no confusion. 

J It should be noted that Jc is positive in the case of air, but it is a complex number in 
the interior of the Earth. 
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The boundary conditions which express the continuity of the tangential 
components of electric force and magnetic force, when the reflector is taken 
to be a perfect conductor, are effectively 

£(n 0 + n*)= pilli, (r = a), 

~ (rllo + rO,) = J- (rlU), (r = a), 

a 

5 - (rrio + rll^) = 0 . (r = c). 

dr 

These conditions yield the following system of equations to determine a ny b n 
and c n :— 

(/3/k) [ctrXn (ha) + {% n (kb) + c n } yfr n (M)] = (fiifh) b n y^ n (ha% 
a r&n (ha) + { (hb) + cf '\jr n / (ha) = b n ^ n f (ha), 

{a n + ir n (hb )} f B ' (he) + c^n (he) = 0 . 

On solving we find that 

a n ~ (hb)'if/'n(ha) y^t n (ha) ‘\jr n (ha) 

U (kb) y/m' (ha) (ha)/(0ih) ^ (ha) 

ir H (kb)^(kc) 0 fnfke) 

&n “ Kn (ha) '\^/' n (ha) Kn (hb) fn (ha) 

Kn Qca) 0hfn {ha)Km g, (hb) f n ' (ha) 

6 .'(*c) 0 yfr n (hb)&(hc) 

where A w =j t n (kct) yfr n (ha) \[r n (ha) 

Kn (lea) (ha)Km ir n ' (ha) 

Kn (kc) 0 yj/'yi (he) 

and hence the value of the Hertzian function just outside the surface of the 
Earth can be reduced to 

"I oo 

U b (a,0) = 2 (2» +1) P, (fi) fn (ha) (hb) (he) 

ruCLO % ^ Q 

-ir n '(hc)Z n (m+** 

It should be noticed that this expression does not reduce to the expression 
given in my former paper when e is made infinite; the reason for, this is 
that the presence of a reflector at infinity would produce convergent waves 
while the Hertzian function of the previous paper contained only divergent 
waves outside the sphere r = b. 

3. How that the general form of the Hertzian function has been obtained, 
we must apply a transformation of the type introduced in my previous paper 
in order to render the function susceptible to numerical computation. 

With this object in view, it is desirable to obtain the corresponding formula 
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for the special case in which the Earth is supposed to be a perfect conductor, 
and then to modify the solution by taking account of the fact that /3 TciftfSJc} 
is not zero. 

4. In the case of a perfect conductor, the Hertzian function is 


n 6 <°) (a, 0) = 


t % /9 \ np / (kV) {kc) (kc) (kb) . 

kab n ir U ^ 1} rnW W ( ]ca ) (&0- (kc) Kn (to) ’ 


and it is to be noted that resonance occurs when any of the denominators 
vanish. As it will be found that resonance effects do not occur in the actual 


physical problem (though they may occur in the ideal problem now under 
consideration), it is convenient to disregard them. 


The series which has just been obtained for IV°) (a, 6) is the sum of the 
residues of 


_2stt ' P{fot-3 (kb)^s ~±(faQ &- 

kab * cos S7r (ka) (kc)~~^ s _ (to) (to)} 

at the poles s = J, f, ; and since it is assumed that there are no¬ 

resonance effects, these are all simple poles. 

How those poles of the function which are not zeros of cos sir are zeros of 

V'W ( ka ) &-*' ( hc ) — (kc) (to); 
and it will appear in § 5 that the only zeros of this function of s lie on the 
cruciform figure,* which consists of the imaginary axis in the s-plane 
together with that part of the real axis which joins the points s = + kc. 

We now consider the integral 

^5 f gP «-*(--/0 {'h-l (kb) (kc) - yfr a ^' (kc) Zs- h (kb)} , 

kdb J cos S7r * (to) ( kc(&0 &-*' (to)} ? 

taken round a contour which consists of a semicircle of large radius, E, with 
centre at the origin, lying on the right of the imaginary axis, together with 
that part of the imaginary axis which joins the points + Ri; on the latter 
part of the contour, Cauchy's principal value has to be taken, in view of the 
poles on the imaginary axis. It is readily verified that the integrand is an 
odd function of s, and so the integral along the imaginary axis is equal to iri 
times the sum of the residues at the poles on it, the rest of the integral 
cancelling. On making R tend to infinity through such values that the semi¬ 
circle never passes through a pole, we find (exactly as in my previous paper) 
that the integral round the semicircle tends to zero. It follows that the sum 
of the residues of the integrand at all the poles on the positive half of the 
real axis plus half the sum of the residues at the poles on the imaginary axis 


* The change which the presence of the reflector produces in the topography of the- 
poles is very remarkable. 
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is zero; the integrand being an odd function of s, the latter term may be 
replaced by the sum of the residues at the poles on the negative half of the 
imaginary axis. 

For brevity, we shall write x, X, y respectively for a, b, c ; and we shall 
denote those zeros of 


fa-h (*)£>-§' (y)~i r s 


! 

1 

S 


(; y ) 



which are under consideration by vi, V 2 , v% y ..., so that either E(^)>0 and 
I ( v ) = 0, or else E (v) = 0 and I (v) < 0 ; and we have 


IT & ( 0 > (a, 6) 


9 


7T 


2 j (~/*) {V'wOO &-t (y) 

kab v cos v7r [d{^r e ^ i , (x)^'(y)—^r s _ i \y)^,_ i '(x)}/ds]. 


5. We shall next prove the theorem which has just been assumed, namely, 
that the number of real values of v is finite, that the number of purely 
imaginary values is infinite, and that there are no other complex values of r. 
It will also be shown that the imaginary values of v are such that the series 
just obtained for II^ 0 ) (a, 0 ) converges very rapidly except when 0 is quite 
small; and that the convergence is uniform with respect to 0 and X when 
0 O < 0 =s 7 r and y > X ^ x, where 0 O is an arbitrary small positive angle. It 
follows from these results that term-by-term differentiation with respect to 0 
is. permissible, and also (by Abel’s theorem on continuity) that x may be 
written in place of X to obtain the effect due to an oscillator placed on the 
surface of the Earth. 

The only inconvenient property of the series is the number of real values 
of v ; there may be as many as 60 or 80 such values. The corresponding 
terms of the series are of an oscillatory character, and in this ideal problem 
numerical computation would be somewhat tedious. 

We now define the function F (s) by the equation 


F 0) - (x) £,_/ (y) -( y) (x) 

= 0^ [{Xy)i {Hs<1> ^ Hs<2) (2/) ~ Hs<1) &> Hs<2> (*)>]■ 


and we proceed to investigate the zeros of F (.$). 

As in my previous paper, we write s = a?cosh^, and we shall also write 
5 = y cosh v ; further, we define S/ 1 ) pc) and S s (2 > pc) to be the (discontinuous) 
solutions of Bessel’s equation, whose asymptotic expansions (of Debye’s 
type) when s and x are large but not nearly equal, are given by the formulae 

S s W(x) — e ^nhu~uoo S hu)-im^^y{i Tr^sin (—m)}, 

S s < 2 > pc) — sin ( - in)}. 
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It follows from the Tables given in my previous paper (p. 93) that 


F (s) = i tt * [( ^)i {S .a> (x) B s m {y) _ S .a> (3/) S t « (a?)}], 
oxoy 


except in narrow bands of the plane, in which the Hankel-Neilsen functions 
of x and y have asymptotic expansions of different types, the function of di 
involving two functions of the type S, while the function of y involves only 
one such function. In these bands the expression just given for F (s) has to 
be modified by the insertion of terms which are negligible in comparison 
with the terms already present, unless s is nearly equal to ±x or ±y. 

Hence the only zeros of F (s) in the whole plane are near the curves on 
which the real part of 


x (sinh u—v, cosh u)—y (sinh v—v cosh v) 

vanishes, or else they are near the points ±x, ±y. 

Regarding x and y as fixed ( s , of course, being a complex variable), we 
write 

y ( 5 ) = x (sinh u—u cosh u)—y (sinh v—v cosh v). 


It is found that %'( s ) = v—% and that R%(s) vanishes on the cruciform 
figure consisting of the imaginary axis and of that part of the real axis 
which joins the points s = ±x. 

Also, when | s ] is large compared with x and y , 

“%(«) = slog(y/x)-i(y 2 -x 2 )/s-«... 

It is now easy to prove that R%(s) does not vanish except on the 
cruciform figure. For if there were any curve on which R^ (, s) vanished, 
this curve would either (i) be closed, or (ii) have points at infinity, or 
(iii) terminate on the cruciform figure, or (iv) terminate on the cuts from 
+ x to + 00 . 

These possibilities do not in fact occur, for (i) the closed curve would 
contain a zero of y? (s), and there is no such point except s = 0 ; (ii) R^ (s) 
vanishes at infinity only on the imaginary axis; (iii) the end points would 
be zeros of y (<§), and no curves on which R% (< 9 ) vanishes emerge from the 
origin except the axes; (iv) lly (s) is negative on the cut from -f x to + co , 
and it is positive on the other cut. 

We shall next determine the number of zeros of F (s) in an arbitrarily 
large region of the plane, and we shall show that this number is the same as 
the number of zeros of F (s), which lie on the cruciform figure formed by the 
imaginary axis and that part of the real axis which joins the points ±y: 
and it is then obvious that all the zeros of F (s) lie on this cruciform figure. 

Consider | {cl log F (s)/ds}cls taken round a rectangle in the s-plane with 
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corners at the points ±M+ (N 4 - J) 7 ri/log when 1ST is a large integer. 
When we make M -► oo , the integrals along the vertical sides of the rectangle 
tend to zero. 

Since N is large, we may replace the functions S 5 <2) by the dominant 
terms in their asymptotic expansions, and so the change of phase of F(s) 
when s traverses a horizontal side of the rectangle is the same as the change 
of phase of sinh x (s) } hence the total change of phase of F (s) as s describes 
the rectangle as 4 Ntt. The total number of zeros inside the rectangle is. 
therefore 2N. 

Next consider the zeros of F (5) which lie on the cruciform figure formed 
by the imaginary axis and the part of the real axis joining ±x; on this 
curve x ( s ) is a pure imaginary. Also, by using Debye's integrals, it may be 
shown* that, with each value of s, which makes tt equal to an integer,, 

is associated one and only one zero of F (5), and this zero is on the cruciform 
figure; the difference between the value of s, which makes ix(s)Jir an 
integer, and the corresponding zero of F (5) may be calculated by Stokes' 
method when u is not small. 

Now ix (#)/ 7 r = {y 2 --x 2 )%-~y Q,o&~ l {xly)}j7r) 

we take Ni to be the greatest integer contained in this number, and write 
N = Ni-f N 2 . Then, as s passes along the real axis from x to 0, and then 
along the imaginary axis to (N-f 7ri/log(yfx), ix( s )/^ increases from 
iX 0*0/^ to (approximately) N +J, and so s passes through N 2 zeros of F (s). 

To show that F (s) has Ni zeros between s = x and s = y, we express F ($) 
in terms of Bessel functions of the first and second kinds by the formula 

m = \ •ni dg- {T* (*) J s (y)—T s (y) J s (*)}]. 

Now as s decreases from y to x, the function defined by the equation 

( r \ — ( i (°°) \/ x J Id {Js (f) \/r} 

y j S — ~ / j - » 

ax i ax 

remains bounded,! since the order of the Bessel function exceeds its 
argument; but © 5 (y) passes through an infinity for each zero of its 
denominator, and there is one such zero for each multiple of 7 ri, which 
3 /(sinh v—v cosh v) passes through; since this function varies from 0 to ix( x ) 
through purely imaginary values, there are Ni such zeros. 

Also © s ( y ) is a decreasing function of s , and it may be shown (by using 
Bessel's equation) that © y (£) is an increasing function of £ when £ < y, so 

* The proof is substantially identical with that given in the simpler case of a single 
Bessel function (of the first or second kind), 4 Boy. Soc. Proc., ? A, vol. 94, p. 193 (1918). 

t ‘Proe. Bond. Math. Soc.’ (2), vol. 16, p. 166 (1917). 
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that % y (x) < % y (y). Consequently the graphs of 0 S (x) and 0 s (y), qua 
functions of s, intersect exactly Hi times when s lies between x and y, and 
each intersection is just on the left of the corresponding value of s which 
makes y (sinh v—v cosh v) a multiple of 7 ri. 

Hence there are Ni zeros of F (s) between s = x and s — y\ and since 
F (s) is an odd function there are evidently 2 (Hi + N 2 ) zeros in all on the 
(large) cruciform figure. 

Hence all the zeros of F (s) are on the cruciform figure; and this is the 
result which had to be established. 

The number of positive zeros of F ( s ) is the integer next less than 
1 % ( 0 )/ 7 r. Since i% (0) = y—x, we can easily construct the following Tables. 
In these Tables h denotes the height of the reflector above the surface of 
the Earth, and X denotes the wave-length, both being measured in kilometres. 

X = 5. 

h 100 150 I 200 


( 0 )/ 7 r 40 60 j 80 


X = 10. 


h , 

100 

150 

200 

%X (0)1-rr 

20 

30 

40 


6 . We now consider the value of 


Wv -1 (X) £,-4' (y) - fr-i 0 v ) 5^4 (X)} /r (v). 


When x, y> and X are regarded as fixed, while v tends to infinity through 
its purely imaginary values, it is easy to show that this function is 
bounded, and so the series for the Hertzian function converges like 
S^sec vtt . P J ,_i (—cos 6) ; and in view of the rapid convergence of this series 
(except near 6 — 0) it follows that the series for the Hertzian function con¬ 
verges uniformly when x^lLsy. 

The statements made at the beginning of § 5 are therefore established. 
We shall consequently suppose that the transmitter is actually on the surface 
of the Earth; the error produced by neglecting the height of the transmitter 
is minute. 

Furthermore, if we use the approximations 


as; 1 ) (a?) 
3s 


-uSP(x), 



3 s; 1 ) G) 

dx 


sinh u . S; 1 ) (x), 


dS s W(x) 


sinh^. S s W(x), 


which are valid when S is not nearly equal to x, we find that 

£,_ 4 ' (I y) - •'h-/ (y) £,-4 (*)} /F' (v) ~ 1 / [(«-«) sinh u] s = „ 
and this is greatest when v—u or sinh u is very small, i.e. } when v is nearly 
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equal to 0 or x. The residue is not really infinite when v is 0 or x, because 
closer approximations have to be used in calculating it. 

Hence n a (°) ( a , 0) is expressed as a series which consists of a finite number 
of terms which are oscillatory functions of 6 combined with an infinite 
number of terms each of which is roughly proportional to an expression of 
the type cosec 6 . exp (~ | v j 6). 

It seems unnecessary to examine this series more closely; our next object 
must be to consider modifications produced in the oscillatory terms by taking 
into account the imperfect conductivity of the Earth. 

7. With the aid of the results contained in §§ 4-6 we shall examine the 
problem in which it is supposed that the Earth consists of an imperfect 
conductor surrounded by a perfect reflector in the upper regions of the 
atmosphere. 

The contour integral to be considered is now 


where 


f sP,_j, (— p-) (hb) £_t' ( lie) — yfr s _ i ' (ha) g,_ j, (kb) d 

hob J COS S7T G($) 


G (s) = (lea) £ s _j' {lcc)—ty s W (he) (ha) 

~m { * s ~i {ka) &-*' W &-* (*“)}■ 

PiE Y' ^ i yh iCt) 


Now in the case of signals over a sheet of water the high conductivity of 
the water ensures that the disturbance at any depth is negligible, and so 
the solution is the same as if the Earth consisted of a sphere of water; and 
in the case of signals over dry land the part of the disturbance which 
does not travel more or less directly along a geodesic from the receiver to 
the transmitter, produces little effect on the transmitter and so the solution 
is the same as if the Earth consisted of a sphere whose conductivity is that 
of dry earth at the surface of the Earth. 


NOW Ip 2 IB — Co) + €{, 

where co is the frequency of the waves and the values of <n, e t -, and cr t */&> are 
given in the following Table :— 



07 . 

< 7 - 

(Til lt>. 

Earth . 

107 

4 

5 *27 A. 
225,000 A 

Water . 

4 *26 x 10 11 

81 



It follows that, in both cases, | Ba | is large compared with the real values 
of r, and the phase of kia is nearly equal to — \ir ; also 


(3h/(/3i!c) — k/ki, 


x 


VOL. XCV.—A. 
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and this is quite small. Hence the zeros of G (s) are associated with zeros of 
F (s), and in calculating the zero of G (s) which is associated with any 
particular real zero of F (s), we may take 

— % ‘yj/'s—l ^ i. 

Let r-f 5 be the zero of G (s) which is associated with the zero v of F (s). 
In estimating the value of 5, we have to use different approximate formulae 
for the various Bessel functions according as v is or is not remote from x 
and y. 

In the easy case when v is not near x we have approximately 


(s -J- 5) — ^ (s) -j- (v u) 5, 

and hence 5 is determined by the approximate equation 

sinh w sinh + u) 5} = (i&/& t *)cosh {% (v) + (v—u) 5}, 
where sinh^(^) = 0. Hence 


5 = ~—, loa sinli u + ikjki 
2(v—u) & sinh u — ikjki * 

Now v—u is a pure imaginary when v lies between 0 and x ; and so 1(5), 
which is our essential requirement, is given by the formula 

sinh u + ikjhi 
sinh u—ikjki 

Except near the critical points 0, x , y , it may be shown that (v—u) sinh u is 
the change produced in cosh u when x is replaced by ?/. Hence 


i(5) 




2 v— 


u 


log 


(v— w)sinh^ = (yjy) — ( v / x ) ^ —(h/a)coa\\n, 
where h is the height of the reflector. It follows that 

‘ 1(5) = 


a 

tanh?/ 

. loo; 

sinh u + ikjki 

2h 


O 

smh u—ikjki 


For simplicity of calculation we shall take the phase of k/h to be exactly 
\ir, and then we shall write 


u = i<jb, k = k{ 8 exp \iri. 


We deduce that 


1(5) = 



tan <f) log 


8 2 -f- Sy/2 sin <f> 4- sin 2 
S 2 — hy/2 sin (j> + sin 2 <j) ’ 


The expression on the left increases steadily from* — oo to 0 as <£ decreases 
from \i r to 0, but its rate of change is, on the whole, quite slow. Thus, so 
long as sin <j> is fairly small (but large compared with 8), the expression is 


* When v — 0, IfS 1 ) is not actually infinite, because our approximation then becomes 
inadequate ; and it turns out that in reality 1(5) is only moderately large. 
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nearly equal to —a$/h^/2); and when sin <£ = S it is nearly equal to 
— 0 ’4:4: ah/h. In the former case we have approximately 


a = 


ikjki 

(v—u) sinh u 


—aikj{hki); 


so the phase of $ is —\nr, approximately. 

In the ease of values of v which are nearly equal to x or y , we have to use 
the appropriate expansions* for Bessel functions of large order. As a 
moderately simple case which can be worked out in detail, we shall consider 
what happens if a value of v is actually equal to x. 

We then have 


F(x) = {H s (1) (x) H s < 2 > (y) —H s (1 > (y) H s (2) (*)}].,= XJ 

= 0 

G (a + 5) = * tt [JL (ay)* {H.W (a) H.W (y)-H s « (y) H s < 2 > (*)}],=,+, 

-frfik/ki) (ay)* {H s <« (a) H s < 2 > (y)-H.« (y) H s ( 2 > (a) 

= 0 . 

Also H, + ,« (y) = ( y ), H I+4 « (y) = e*»H,«(y), 




where ^ has the value corresponding to s = a?, so that v = if{2hja). 

On using these results to simplify the equation G(^ + d) = 0, and then 
eliminating the Bessel functions involving y by means of the equation F (x) = 0, 
and making the further simplification produced by the approximate equation 


a 


^ {xi H s (1 > (*)} 


_ e %m 


0 


~ (a)} 


!s=x 


(which follows at once from a formula due to Debye), we find that 

~ [xi {FL S {1) (®) (#) e -1 *} ]«=*+$ 

-(ik/h) [xi {H s ^ (a?) (*.) = 0. 


If we put | = %x~i( — 23)* and substitute for the Bessel functions their 
approximate values in terms of Bessel functions of orders -h and —this 
equation may be written 

(3£/a)* [BL** 1 ) -*<»>(£) 

- (iA/A,-) [EL^ (£) (f) = 0; 

that is to say, 

(3 f;/^)* [J_|(^) sinh ?;3 -f J| (f) sinh (vS- -4- ± 7 n)] 

-{•{ikjki) sinh (vS- + ^7ri )—Ji (£) sinh ^3] = 0. 

* 4 Proc. Camb. Phil. Soc., 5 vol. 19, p. 110 (1917); ‘Proc. Boy. Soc., 5 A, vol. 95, p. 96 
(1918). 
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Now if the value of 3 given by this equation is not appreciably larger* 
than the value previously found, £ does not exceed about \ (when the wave¬ 
length is short), and it is less with a greater wave-length. Hence, in approxi¬ 
mating to the Bessel functions, we may safely neglect all except the first terms 
of the functions of negative order, so that 


sinh v$ 

sinh (^3 + ^ 7 ^) 


(ik/ld). Qx)i r (j)/r (f) = 7), 


say; and, to a sufficient approximation, 

v3 = ^ v / 3/( 2 — rj), 

where v = i^/(2hja) and \2—rj\<2 } 

since arg rj = — r. 

Hence the modulus of $ is greater than if 5 were given by the formula 

vS - i Vs /( f), 


and the phase of 3 is less than — 

Now the imaginary part of 5 given by this formula is numerically greater 
than that given by the formula obtained when v is not nearly equal to x if 

V(! a/ A) • (I *) 4 r a)/r(|) > 0-63 (a/h), 
i.e ., if h > 4*9 AJ. 


Now the greatest value of wave-lengths used in practice is 10 kilom., and, 
since it is certain that h> 23, the zeros of G(s), which are nearly equal to x, 
are less important than those whose real part is appreciably less than x. 

The zeros of G ( s ), which are nearly equal to y, may be treated in a similar 
manner, with the same result. 

Now, in the case of waves over dry land, 8 is about as small a value as 
may be given to sin </> when the simple approximation is used; in the case of 
waves over water, the smallest admissible value of sin </> is much greater 
than S. 

Hence we may infer that, in the case of waves over dry land, the bulk of 
the zeros of G (5), which are in the neighbourhood of the real axis, are 
between the lines on which I($) has the values —0*7 a8/h and 0*44aS/A, 
none being nearer than the latter line ; in the case of waves over the sea, 
they are clustered much more closely about the former line, except the zeros 
near the point x ; these zeros are nearer the real axis when the wave-lengths 
are long (but they are not nearer in the case of the comparatively short 
wave-leugth of 3*75 kilom. used by Austin; for a wave-length of 10 kilom., 
they are not nearer if h > 60, as is probably the case). 


* If the value were appreciably larger, the contribution to the Hertzian function of 
the corresponding term would be smaller. 
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The values of 3 corresponding to different values of v are indicated roughly 
by the figure. 



It should be observed that no zero of Gr ( s ) can lie on the positive half of 
the real axis; for, if it could, then either (i) kc 1 ^-1 ' (ka)/f s - i (kid) would 
be real with 5 positive and hi complex, and this is not the case ; or else (ii) 
both 

fa-h (°°) (y) 

and its derivate with respect to % would vanish simultaneously; this would 
involve %s-h (y) and tys-i (y)> having a common zero, which is not the case, 
or else the Wronskian of yjr s ^(x) and £ s -.$(x) would vanish, which is impos¬ 
sible, since its value is —i. 

8 . If, now, we take the integral introduced at the beginning of § 7, and 
evaluate it in the usual manner, we find that the transformed Hertzian 
function consists of three parts. The first part is a finite series of not more 
than 80 terms, corresponding to the real zeros of F ($); in this series the 
terms are of the order of magnitude of P„_a (— cos 9) sec. vnr, where I {v) 
lies between — 0 ,Y 1 aSjh and —0*44 aS/h; these terms may be taken as 
equivalent to y'eosec 9 . exp ( — ^aZOjh) in absolute magnitude, except when 
9 is nearly equal to ir. 
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The second part is a rapidly convergent infinite series, in which the terms 
are of the same type as in the first series, but v tends to infinity in the 
direction of the imaginary axis, so that the Legendre functions may be 
replaced by negative exponentials. This series is negligible compared with 
the first series.* 

The third part is the integral 


2tt p* sPss(-fi) {yfr s ^j.(1cb) ga-V (kc) — jr s __i (he) (kb)} 


hob J, 


• oo % 


cos sir 


G(«) 


ds 


2 tt 


cos sir 


{■ .* m w Qcc) w wi 


Now 1/G (s) — 1/G ( — s) contains as a factor 


1 


1 


.G (s) G (—s)_ 




and, as in my former paper, the presence of this factor makes the integral 
(or its principal value) negligible in comparison with the first of the two 
series. 

Hence the order of magnitude of the Hertzian function is 

y/ cosec 6 . exp ( — | aS 6 j h), 
except near the antipodes of the oscillator. f 

Now \ahjh = \a(<Tij co)~^/h 

_ 10 4 / 27 rC\^ 
irh \ (Ti\ j 

In the case of waves over dry land, this is 1400 /(hy/\); in the case of 
waves over the sea, it is 7 j(h^/X). 

The latter expression is evidently excessively small compared with the 
value suggested by Austin, 9'6/vA. The inference is that the upper regions 
of the atmosphere cannot be treated as a perfect conductor. 

We shall therefore have to attack the problem of an imperfectly 
conducting sphere surrounded by a concentric shell of imperfectly con¬ 
ducting material, the region between the sphere and the shell being filled 
with homogeneous dielectric.! 


* There is no sharp dividing line between the end of the first series and the beginning 
of the second. 

f Near the antipodes the factor >v / cosec 6 has to be suppressed. 

t The investigation of the terms of the Hertzian function corresponding to zeros of 
G($) near s = 0, and between s = oo and s = y, has been suppressed deliberately, in order 
that the paper should not be overloaded with analysis. 
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9. In the problem which has now to be considered, we suppose that the 
constants in the upper atmosphere have the values /3 e , y e , /c e ; since the upper 
atmosphere is a good conductor, the Hertzian function is inappreciable when 
r exceeds c by a comparatively small amount, so that the error produced by 
supposing that the conductor extends to infinity is negligible. 

We take 11^ and Hi to have the same forms as in § 2, and we suppose that 
the Hertzian function in the conductor is H e where 

• oo 

Ue ~ ( 2n +1 ' ) dn ^ n ^ p ” (jx) ' 

The boundary conditions are now 

/3(IIo + n d ) = (r = a), 

a - (rllo + Tll d ) = £ (rlli), (r = a), 
or or 

^(Ilo + IIrf) = (r = c), 

(rn„ + rU d ) = ^ ( rUe ), (r = e). 

The equations connecting a n , b n , c n , d n are 

(^ jk) \_^nKn (Jed) ~b (Jeb) -f- C n ) (led )J — (j3 jlej) ( h{Ct\ 

Q"rX,n (Jed) ~ J- { (Jeb) T” ^ri\ (Jed) b n / \jr n (JeiCt) i 

08/*) (&<^)} K>n (^0 T* Cntyn (^)J =: ifejl^e) $v&n 

T* 'ty'n (kb )} T" (kd) == (kgC). 

On solving these equations we find that the value of the Hertzian function 
at the surface of the Earth is 

leab ^ ^ (/ l )/h n , 

g n = ^ (kb) £n (hc) — y\r n ' (Jec) (*&) 

-f| {*. W 5. (to)-*„ (to) S. (to)}, 

hn == (led) (kc) (kd) (Jed) 

-f| W”"' <*“> &• ( fo ) &' ( te ) 

+ || Hn (to) & (to)-*. (k*) 5. (to)}]. 


where 
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The value of the Hertzian function is consequently dependent on the zeros 
of h s ^. If we denote the zero of h s _% which corresponds to the zero v of F(s) 
by v + $ (so that $ has not the significance as in §§ 7-8), the equation which 
determines $ approximately is 

sinh u sinh v sinh (v) + (y—u) —(ilc/ki) sinh v cosh (v) + (v —u) 

- (ik/k e ) (l - [sinh u cosh {% (v) + (v-u) 3} 


— (ikjki) [sinh {%(^) + G— v) 3}] = 0, 

where sinh %( v ) ” 0 and s = v + §. On solving we find that 
(v—tc) § = tanh -1 [(ik/ki) cosech u] 


q-tanli 1 [(ilc/lc e ) (1 




cosech v\ 


Hence $ consists of two terms each of which is of much the same type as 
the single term investigated in § 8. 

In order to obtain an expression of the form taken empirically by Austin, 
it is necessary that h e should be large compared with h, and then the phase of 
lc e is nearly — r; and in the case of waves over the sea we obtain Austin's 
factor 9*6 if 

0*44 (a/h) . \ (k/Jc e )\ = 9*5 V*, 
i.e., if 0*44 (a/h) ( 27rC/<j e )i = 9*5. 

This equation gives h?cr e = 1*67 x 10 11 . 

Hence Austin's factor is obtained if the reflector has a height of say 
100 kilom. and conductivity 1*67 x 10 7 . A conductor with this conductivity* 
is a fairly good conductor so that Austin’s result seems to be quite 
compatible with what is known concerning the upper atmosphere. The 
factor corresponding to Austin’s number 9*6 in the case of waves over dry 
land is 9*6 + (1400/A). I am not aware of any quantitative experiments 
having been carried out for such waves, but the increase in A is quite 
consistent with the well-known fact that electric waves are transmitted 
much more easily over water than over dry land. 

The fact that h s _ i has a number of zeros with the same imaginary part, so 
that the Hertzian function contains a number of terms of the same order of 
magnitude but of different phases, suggest that, in certain circumstances, 
interference effects may be produced; such effects have actually been observed. 


* The measure of this conductivity in ordinary electro-magnetic units is l'44x 10~ 15 ; 
the corresponding specific resistance is 6 95 x 10 14 . For copper the corresponding 
resistance is about 1600 ; for fresh, water about 2 x 10 14 . 
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10. It is to be concluded that the ionisation theory which results in the 
upper regions of the atmosphere acting as a reflector of electric waves of long 
wave-length is sufficient to explain the observed facts concerning the rate of 
decay of electric waves transmitted over large sheets of water; and it has 
been seen in my previous paper that a theory of pure diffraction is insufficient 
to account for these facts. The results contained in this paper consequently 
afford a confirmation of the theory put forward by Heaviside and others 
and modified by Eccles. 
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